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generated mass to the one obtained by means of an average site occupation constraint. 
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I. INTRODUCTION 

Quantum Electrodynamics QED(^2+i) is a common 
framework aimed to describe strongly correlated systems 
such as quantum spin systems in 1 time and 2 space di- 
mensions, as well as related specific phenomena like high- 
Tc superconductivity 9, 10, 20]. A gauge field formula- 
tion of antiferromagnetic Heisenberg models in d = 2 
dimensions leads to a QED3 action for spinous, see f. i. 
Ghaemi and Scnthil 9] and Morinari lOj. This descrip- 
tion raises the problem of the mean-field solution and 
the correlated question of the confinement of test charges 
which leads to the impossibility to determine the quan- 
tum fluctuation contributions through a loop expansion 
in this approach [2ll2ll23l|. 

We consider here the 7r-flux state approach introduced 
by Affleck and Marston |^ 0] . The occupation of sites 
of the system by a single particle is generally introduced 
by means of a Lagrange multiplier procedure 0, Q ■ In 
the present work we implement a strict site-occupation. 
It can be constructed by means of constraints imposed 
through a specific projection operator which introduces 
an imaginary chemical potential. This has been proposed 
by Popov and Fedotov [l| for SU{2) spins and generalized 
by Kiselev et al. ^3 to SU (N) semi-fermionic Hamilto- 
nians. It is our aim in the present work to confront the 
outcome of the two approaches. 

Here we concentrate on the behaviour of the spinon 
mass which is generated dynamically by an U{1) gauge 
field. Appelquist et al. 0, showed that at zero tem- 
perature the originally massless fermion can acquire a dy- 
namically generated mass when the number N of fermion 
flavors is lower than the critical value Nc = 32/7r^. Later 
Maris confirmed the existence of a critical value 
Nc — 3.3 below which the dynamical mass can be gen- 
erated. Since we consider only spin- 1/2 systems, N = 2 
and hence N < Nc- 
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At finite temperature Dorey and Mavromatos [13 and 
Lee 01 showed that the dynamically generated mass 
vanishes at a temperature T larger than the critical one 
Tc. 

We shall show below that the imaginary chemical po- 
tential introduced by Popov and Fedotov [J modifies no- 
ticeably the effective potential between two charged par- 
ticles and doubles the dynamical mass transition tem- 
perature, in agreement with former work at the same 
mean- field level ,17]. 

The outline of the paper is the following. In section UTI 
we recall the projection procedure introduced by Popov 
and Fedotov (PFP) leading to a rigorous constraint on 
the lattice site occupation. In section HTll we derive the 
Lagrangian which couples a spinon field to a U{1) gauge 
field. Section llVl is devoted to the comparison of the effec- 
tive potential constructed with and without strict occu- 
pation constraint. In section^Jwe present the calculation 
of the mass term using the Schwinger-Dyson equation of 
the spinon. 

II. SITE OCCUPATION CONSTRAINT FOR 
QUANTUM SPIN SYSTEMS AT FINITE 
TEMPERATURE. 

Heisenberg quantum spin Hamiltonians of the type 

H=-^JijSiSj (1) 

with {Jij} > can be projected onto Fock space by 
means of the transformation 

St = flkl 

Si = ■^{fli;fi,T " fliAi) (2) 

where {fi,a,fi,cr} are anticommuting fermion operators 
which create and annihilate spinon with cr = ±1/2. 
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This transformation is not bijective because the di- 
mensionahty of Fock space is larger than the dimension- 
ahty of the space in which the spin operators {5*^} are 
acting. Indeed, in Fock space, each site i can be occu- 
pied by 0, 1 or 2 fermions corresponding to the states 
|0,0 >, |1,0 >, |0, 1 >, |1, 1 > where |0,0 > is the parti- 
cle vacuum, |1, >= | + 1/2 >, |0, 1 >= | - 1/2 > and 
|1, 1 >= I -I- 1/2, —1/2 > in terms of spin 1/2 projections. 
Since one wants to keep states with one fermion per site 
the states |0, > and |1, 1 > have to be eliminated. This 
can be performed on the partition function for a system 
at inverse temperature /3 

Z = Tr [e-^"] 

where the trace is taken over the whole Fock space by 
the introduction of a projection operator 



where 77 is a lattice vector {01,02} in the Ox and Oy 
directions the Hamiltonian takes the form 



<ij> 

where i and j are nearest neighbour sites. 

The number operator products {niUj} in Eq.(|3J) are 
quartic in terms of creation and annihilation operators 
in Fock space. In principle the formal treatment of these 
terms requires the introduction of a mean field proce- 
dure. One can however show that the presence of this 
term has no influence on the results obtained from the 
partition function. As a consequence we leave it out from 
the beginning as well as the contribution corresponding 
to the {rii} terms. 



Tr 



where N is the particle number operator and /i = i7r/2/? 
an imaginary chemical potential [ij . Indeed, the presence 
of the states |0,0 >i and |1, 1 >i on site i leads in Z to 
phase contributions which eliminate each other 



= 



A. Exact occupation procedure 

Starting with the Hamiltonian 

J 



(4) 



the partition function Z can be written in the form 



and hence the contributions of these spurious states are 
cancelled as a whole. 

The common alternative approximate projection pro- 
cedure would be to introduce a chemical potential in 
terms of real Lagrange multipliers {A^} 



Z = Tr 



where Ui is the particle number operator on site i and 
the {Xi} are fixed by means of a saddle point procedure. 



-^({S*,„,C..<.}) 



where the ^1,0-} are Grassmann variables corre- 

sponding to the operators {fj„, fia} defined above. They 
depend on the imaginary time r in the interval [0,/?]. In 
the continuum limit the action A is given by 



•^0 V ^,a 



III. SPIN STATE MEAN-FIELD ANSATZ IN 2D 



+ H({?*,(r),e.,.(r)}) 



In 2d space the Heisenberg Hamiltonian given by 
Eq.(^ can be written in terms of composite non-local 
operators {T>ij} (" diffusons" ) 0] defined as 



If the coupling strengths are fixed as 



where 



n{T) = H{t) - hN{t) 



(5) 



and N{t) is the total particle number operator. A 
Hubbard-Stratonovich transformation on the corre- 
sponding functional integral partition function in which 
the action contains the occupation number operator as 
seen in Eq.ljSjl eliminates the quartic contributions gen- 
erated by Eq.((2l and introduces the mean fields {Ay}. 
The Hamiltonian takes then the form 
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(6) 



The fields {A^j} and their complex conjugates Kij can 
be decomposed into a mean-field contribution and a fluc- 
tuation term 



A,: 



5A, 



The field A^'^ can be chosen as a complex quantity 

The phase is fixed in the following way. Consider 
a square plaquette □ = {i,i + ex,i + ey,i + ex + ey) where 
and Cy are the unit vectors along the directions Ox and 

Oy starting from site i on the lattice. On this plaquette 
we define 




FIG. 1: The two-dimensional spin Brillouin Zone (black area) 
and the lattice Brillouin Zone (whole square). 



H 



mf 



\J\ 



keSBZ 



H 



(8) 



n ' 



which is taken to be constant. If the gauge phase (j)^^ 
fluctuates in such a way that (p stays constant the average 
of A™'^ will be equal to zero in agreement with Elitzur's 
theorem In order to guarantee the SU{2) invariance 
of the mean-field Hamiltonian along the plaquette we fol- 
low 11 HHHIi and introduce 



e*- 4 ( 1) , if rj — fi + Cx 
g-^-f (-1)' if = ^ e 



with 







H 





-fi + A cos jzjk^^ky -iA sin f ^7fcx,fc„+7r 
HAsin^zjk^^ky+n -A* - Acos f Z7fc^,fc^ 



where A = |A"-^|. The Spin Brillouin Zone (SBZ) covers 
half of the Brillouin Zone (see figure^. The 7j:'s are 
defined by 



7fe 



z 



(cos kxtti + cos kya2) 



V 



where Cx and Cy join the site i to its nearest neighbours j. 
Then the total flux through the fundamental plaquette 
is such that = tt which guarantees that the SU{2) 
symmetry of the plaquette is respected [s^l • 

At the mean-field level the partition function reads 



B. The TT-fiux Dirac action 

As already shown in earlier work by Ghaemi and 
Senthil and Morinari the spin liquid Hamiltonian 
© for systems at low energy can be described by four- 
component Dirac spinous in the continuum limit. The 
Dirac action of spin liquid in (2-f 1) dimensions is derived 
in appendix 1X1 In Euclidean space this action reads 



where 



n 



"^f - \J\ 2^ iJ iJ 

' ' <ij> 

+ E l^y-^v + A:f 4 



- (7) 



as read immediatly from Eq. JBJ 

After a Fourier transformation the Hamiltonian iQ takes 
the form 



Se = 



dT 



J° {dr - /i) + Aj'^dk 



(9) 



where A = 2 A cos j is the "light velocity", and {7^^} are 
the Dirac gamma matrices in (2-1-1) dimensions. Spinous 
move in a "gravitational" field and the metric can be han- 
dled in a Minkowskian (or Euclidean) metric assum- 
ing A = 1 without altering the physics of the problem. 
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Since the Heisenberg Hamiltonian is gauge invari- 



ant in the transformation -0 ~ 
has to be written in the form 



the Dirac action 



Se = 



(fr 



(1 - XjO^d-")] a. 



^-ipra [7m (^M - ^Sfl/x)] -Ipfa \ (10) 



FIG. 2: The dressed photon propagator. Wavy lines corre- 
spond to the photon and solid loops to the fermion insertions 



generate an orthogonal tensor basis transversal to the 
photon momentum q^^ 



Here g is the coupling strength between the gauge field 
and the Dirac spinous tp. In the first term cor- 
responds to the "Maxwell" term —jffivf^" of the gauge 



field 



where = 9^ 



d^Gfj^, A is the parameter 



of the Faddeev-Popov gauge fixing term —A [d^^a^) |23| 
and 5^" the Kronecker 5. □ = 9^ -I- is the Lapla- 
cian in Euclidean space-time. This form of the action 
originates from a shift of the imaginary time derivation 
dr ^ dr + fJ, and leads to a new definition of the Mat- 
subara frequencies only for the fermion fields ■0 Tl which 
read then 



UJF,n = UJF,n - fl/l = —(n+ 1/4) 

This modification will induce substantial consequences 
as it will be shown in the following. 



IV. THE "PHOTON" PROPAGATOR AT 
FINITE TEMPERATURE 



Integrating over the fermion fields ip leads to a pure 
gauge Lagrangian Ca — ^ap^^i^"'^ where A^^ is the 
dressed photon propagator from which we shall extract 
an effective interaction potential V{R) between two test 
particles and extract a dynamically generated fermion 
mass. 

The finite-temperature photon propagator in Eu- 
clidean space (imaginary time formulation) verifies the 
Dyson equation 



n, 



(11) 



The detailed calculation of the polarisation function II^i, 
is given in appendix IbI 

Since the system is at finite temperature and "relativis- 
tic" covariance should be kept the polarisation function 
may be put in the form 



where 11^ and Hb are related to 11^ by 11^ = IIi 4 
112 and Hb = Ha. The expressions of Hi, 112 and Ilg 
are explicitly worked out in appendix^ and B^i, 



A 



t- - 

r 



with ry^^ ^ 5f,i, - , 
- {q.u)u^. Here = 



Uf, ^ ^^Qf, and = 

(1,0,0) is the three-vector of 



the thermal bath. 

The dressed photon propagator A^j, is obtained by the 
summation of the geometric series shown in figure |21 and 
reads 



Hi + n2 q^+ Ha 



^^•^ (1 - 1/A)>'^" 



'(g2)2 



(12) 



A. Effective potential between test particles 

The effective static potential between two test particles 
of opposite chages g at distance R is given by 



V{R) 



-g^ I dTAoo{T,R) 



1 



d^q 
27r J (27r)2 



Aoo(g" =0,g)e 



iq.R 



g^ f°° 1 

— / dqqJoiqR).— ~— 

27rJo a^ + Ha 



g2 + n3(m = 0) 

where Jo{qR) is the zero order Bessel function. 
The polarisation contribution lisiq'^ = 0,q) is 
equal to dxlog2 ^coshPq^/x{^^~x)j when tak- 

ing the PFP imaginary chemical potential into ac- 
count. This has to be compared to the expression 

1^ /g^ da; log 2 ^cosh ^q^J x{l — x)^ when the Lagrange 
multiplier method for which A = is used jl2j. 

For small momentum q ^ 0, H3(m = 0) can be identi- 
fied as a mass term {M^^{(3) f. For R > (M,f ^) 
effective potential reads 



the 
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T = 0.; 




R 



FIG. 3: Effective static potential with (full line) and without 
(dashed line) the Popov-Fedotov imaginary chemical potential 
for the temperature T = {0.001, 0.05, 0.5}. 



G-\k) = G(") \k) 



7pG(p)A^,^(fc -p)r^ 



(13) 



where p = {po = u)F,n,P), G is the spinon propagator, 
Ti, the spinon-" photon" vertex which will be approxi- 
mated here by its bare value g^n and A^i, is the dressed 
photon propagator (|12|l . The second term in p3(l is the 

fermion self-energy S, (G-i = G(0) - S). Performing 
the trace over the 7 matrices in equation (|13|1 leads to a 
self-consistent equation for the self-energy 



\ 



\ 



[k-p) 



p2 -I- S(p)2 



(14) 



In the low energy and momentum limit m(/3) — S(fc) 
S(0) the equation H14|l simplifies to 



0.005 0.01 . 015 0.02 0.025 0.03 0.035 A 



FIG. 4: Temperature dependence of the dynamical mass gen- 
erated with (full line) and without (dashed line) the use of 
the Popov-Fedotov procedure. 



(27r 



-A^^(-p), 



1 



p2 -f m{(3Y 



(15) 



If the main contribution comes from the longitudinal part 
Aoo(0, —P) of the photon propagator l(T3|l goes over to 



1 = 



(27r)2 



V{R,p) 



dq- 



qJoiqR) 

'2 + MP^' 



N\I BttRMPP 



P2 + (m = 0) „ + P2 + m(/5)2 



(16) 



Performing the summation over the fermion Matsubara 
frequencies uJF,n the self-consistent equation takes the 
form 



where N = 2 since we consider only S* = 1/2 spins. 

Figure 121 shows the effective potential between two op- 
posite test charges at distance R ^ {MqP)~^. The 
screening effect is smaller when the imaginary chemi- 
cal potential /i is implemented rather than the Lagrange 
multiplier A. By inspection one sees that {MqP)~^ = 
V2(Mo^=0)-i. 



V. DYNAMICAL MASS GENERATION 

We show now how the PFP doubles the "chiral" 
restoring transition temperature of the dynamical mass 
generation. The Schwinger-Dyson equation for the 
spinon propagator at finite temperature reads 



4ttN 



dPP 

(2^)2 



Ptanh/3\/p2 + to(/3)2 



p2 + n3(m = 0) JP^ + mij3) 



(17) 



Eg. ljlTII can be solved numerically with a cutoff A fixed 
at 00 in an anlytical calculation. By inspection of equa- 
tion p7(l and the corresponding result obtained by Dorey 
and Mavromatos and Lee 13] one sees that the imag- 
inary chemical potential used which fixes rigorously one 
spin per lattice site of the original Hamiltonian dou- 
bles the transition temperature. This result is coherent 
with the results obtained elsewhere ^3 where spinous 
are massless. 
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Since the mass can be identified with a superconduct- 
ing gap one can evaluate the parameter r = ^r|^ where 
to(0) is the mass at zero temperature and Tc the tran- 
sition temperature for which the mass becomes zero. 
Dorey and Mavromatos obtained r ~ 10 and Lee 
[T3I I computed the mass by taking into account the fre- 
quency dependence and obtained r ~ 6. We have shown 
above that the imaginary chemical potential doubles the 
transition temperature so that the parameter r is ~ 4.8 
for a/ A = 00 to compare with the result of Dorey and 
Mavromatos and r ~ 3 to compare with Lee's result. Re- 
call that the BCS parameter r is roughly equal to 3.5 and 
the YBaCuO parameter r ~ 8 as given by the experi- 
ment |18||. 
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APPENDIX A: DERIVATION OF THE 
EUCLIDEAN QED ACTION IN (2+1) 
DIMENSIONS 

At low energy near the two independent points k — 
(±7r,7r) -I- fc of the Spin Brillouin Zone (see figure the 
Hamiltonian ^ can be rewritten in the form 



VI. CONCLUSION 

We mapped a Heisenberg 2d Hamiltonian describing an 
antiferromagnetic quantum spin system into a QED(^2+i) 
Lagrangian coupling a Dirac spinon field with a U{1) 
gauge field. In this framework we showed that the im- 
plementation of the constraint which fixes rigorously the 
site occupation in a quantum spin system described by a 
2d Heisenberg model leads to a substantial quantitative 
modification of the transition temperature at which the 
dynamically generated mass vanishes in the QED(^2+i) 
description. It modifies consequently the effective static 
potential which acts between two test particles of oppo- 
site charges. 

The imaginary chemical potential 0| reduces the 
screening of this static potential between test fermions 
when compared to the potential obtained from standard 
QED(^2+i) calculations by Dorey and Mavromatos 
who implicitly used a Lagrange multiplier procedure in 
order to fix the number of particles per lattice site 0, E| 
since A = at the mean-field level. 

We showed that the transition temperature to "chi- 
ral" symmetry restoration corresponding to the vanish- 
ing of the spinon mass m(/3) is doubled by the introduc- 
tion of the Popov-Fedotov imaginary chemical potential. 
The trend is consistent with earlier results concerning 
the value of Tc 0|. It reduces sizably the parameter 



determined by Dorey and Mavromatos |l^ 



2m(0) 

and Lee 

Marston [s^l showed that in order to remove "forbid- 
den" U{1) gauge configuration of the antiferromagnet 
Heisenberg model a Chern-Simons term should be natu- 
rally included in the QED3 action and fix the total flux 
through a plaquette. When the magnetic flux through a 
plaquette is fixed the system becomes 27r-invariant in the 
gauge field and instantons appear in the system. This 
is the case when the present non-compact formulation of 
QED3 is replaced by its correct compact version [32j . 

It is our next aim to implement a Chern-Simons term 
|3l| in a system constrained by a rigorous site ocupation. 



H 



T3 



keSBZ <^ 

-/il + \/2A 



y2A 



no T3 

T2 

-T2 



kylL 



ky.iTL 




with 7? = (tTjTt) the Brillouin vector. ti,T2 and T3 are 
Pauli matrices 



1 

1 



T2 





i 



, T3 



1 

-1 



f ^ - and f, r ( f ^ - and f„ r ) are fermion cre- 
ation and annihilation operators near the point (tt, tt) 
((-tt.tt)). 

Rotating the operators 




leads to 



keSBZ 



Ak4 



Tl 



Ak^ 



T2 



"+ » T2 J V ^1 

where k^ — k^ + ky and fc_ — k^ — ky, A = 2 A cos ^ and 



^ka 



( f la.ka- \ 
f lb,ka 
f 2aka 
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In the Euclidean metric the action reads 



keSBZ 
T2 



{dr - 



T3 
T3 



Tl 



^ > -Tl / \ -T2 



Through the unitary transformation 



1 W 1 

e'T^3 ) • 1 -Tl 



and writing = fc2 and fc_ = fci 



Se ^ dr J2 H'^fc^ 

keSBZ <^ 

where '0 — and the 7 matrices are defined as 



T3 
-T3 / ' 



Tl 

-Tl ; ' 



T2 

-T2 



Using the inverse Fourier transform -0^^. — J (Pripr< 
the Euchdean action finally reads 



ik.r 



ijja (fcl 



igj^a^iki - ^2) 



(27r)2/3)2 



i>a (^2) 



with k = (oji? = ^{n + 1/4), A:). Integrating over the 
fermion field tJj and keeping the second order in the gauge 
field leads to the effective gauge action 



^eff [a] = -Tr [GF-igj^'a^] 

with Tr ^^-E / $^t'>^- The trace tr extends 

over the 7 matrix space, and Gj^"'^(fci — fc2) — i {2j)2p ^{^i~ 
/c2)- The pure gauge action comes as 



elf [a\ = -g- 



, 1 [ d% 1^ f (Pit' 



2/3^^7 (27r)2-/3^y (27r)2 



^2 .Yap,{ki-k ). ^ .^"ay (-(fci - k ) 



.2 



With the change of variables ki — k = q and ki = k 



Se^ dr d^rJ2^p- k - ^) + 



2/3^7 (2^)2 



a^(-g)n^''(g)a,(g). 



With a "light velocity" = (1, A, A). The covariant 



derivative which takes into account llll reads 



where q = {ujb — ^"^j o) a-nd the polarisation function 
is given by 



where LL>a,ab — {daGvb — T^^e^f,) , are the vierbein 
[2^ for which the metric is defined as (7'^" = ri"^"'e!^e!^ = 
w'^J^'^ with rf" =_-l,^*^ = (5*^ and T is the Christoffel 
symbols. Since A is constant we see clearly that the 
vierbein are also constant, tOa^ab = in a dilated flat 
space-time with the Euclidean metric 17^^ — v^6fj_^. 



APPENDIX B: DERIVATION OF THE PHOTON 
POLARISATION FUNCTION AT FINITE 
TEMPERATURE 

The Fourier transformation of the spinon action given 
by Ea. lfTUI) reads 



fc2 ■ ' ■ ' [k + q) '\ 



Then using the Feynmann identity \ 



/o '^^ (ax + (l- 



(ax + (l-x)by 



n'^" can be rewritten as 



dx- 



kpikfj Qri) 



[{k + q)^x + {I - x)k^Y 



By means of a change of variables k 
k — xq and using the identity tr [7''7'^7''7''] 
4. [Sp^.Srii, - 5pn-5py + Spy.S^j^] one obtains 
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2k ^kj^ 



+ {'^-'2■x){k^q^, + q^,k^)-x{l~x)2q^q^ 



n' 



00 



a I' I <Pk 
— I ax 



(1 ~ 2x)qoS' 



Si -2 



k' + x{l — x)ql 



^fc^ + (1 ~ '2x)k^qfj — x(l — x)q^^^ for the temporal component and 



'2 ol^ 

/ k + a;(l - x)q^ 



N=2 



where a = Following Dorey and Mavromatos 

[12, Lee [13, Aitchison et al. ^ and Gradshteyn [l^ 
we define 



(1 - 2x)qoStj.S 



2x{l - x){q^5ij - qiqj)S2 



for the spatial components. 

Integrating over the fcrmion momentum k' one gets 



^1 - E 



5"* 



ii=-oo ^k'^ + x{l - x)q^ 
\ sinh(27rr) 



E 



cosh(27ry) - cos(27rX)_ 

1 _ /32 1 dSi 



-oo ^k"^ + x{l ~ x)q'^ 

UJp 



r ^2 

n=-oo A;' + a;(l — x)(7^ 



'47r 9X 



with X = x.m+1/4 and Y = k' + x{\ - x)q'^ . The 
polarisation can be expressed in terms of these sums and 
reads 



nO" - na - 4ni - 



where 



aq 



q- 

Hi ( <5y - ^ ) + n2(5y 



da;\/ x(l — x) 



Hi 

n2 



Ha = -r / dx\og2D{X,Y) 



sinh /3(7y^ x{1 — x) 



DiX,Y) 
cos 2nxm 



cum / , , „ , 



and = cosh (/3gy a;(l — x)^ + sin(27rxm). 
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